Parametric derivation of the observable relativistic periastron advance for binary 

pulsars 
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We compute the dimensionless relativistic periastron advance parameter k, which is measurable 
from the timing of relativistic binary pulsars. We employ for the computation the recently derived 
Keplerian-type parametric solution to the post-Newtonian (PN) accurate conservative dynamics 
of spinning compact binaries moving in eccentric orbits. The parametric solution and hence the 
parameter k are applicable for the cases of simple precession, namely, case (i), the binary consists of 
equal mass compact objects, having two arbitrary spins, and case (ii), the binary consists of compact 
objects of arbitrary mass, where only one of them is spinning with an arbitrary spin. Our expression, 
for the cases considered, is in agreement with a more general formula for the 2PN accurate k, relevant 
for the relativistic double pulsar PSR J0737-3039, derived by Damour and Schafer many years ago, 
using a different procedure. 
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I. INTRODUCTION 

A detailed knowledge about the dynamics of spinning 
compact binaries is of interest to both astrophysicists and 
theoretical physicists. The astrophysical interest arises 
from the highly accurate radio observations of relativis- 
tic binary pulsars, which can be used to test general rel- 
ativity in strong field regimes [1]. Moreover, inspiralling 
spinning compact binaries are promising sources of grav- 
itational radiation for the ground based and proposed 
space based laser interferometric gravitational-wave de- 
tectors [2]. The theoretical interest is due to the possibil- 
ity of constructing numerically spinning compact binary 
configurations in full general relativity [see Ref. [3] and 
references therein]. Further, it is of some interest to ex- 
plore if the dynamics of spinning compact binaries will 
be deterministic or not [4, 5]. 

In many of the above mentioned scenarios, the dynam- 
ics of spinning compact binaries can be described in great 
detail by the post-Newtonian (PN) approximation to gen- 
eral relativity. The PN approximation allows one to ex- 
press the equations of motion of a compact binary as 
corrections to Newtonian equations of motion in powers 
of (v/c) 2 ~ GM/(c 2 R), where v, M, and R are the char- 
acteristic orbital velocity, the total mass and the typical 
orbital separation of the binary. For spinning compact 
binaries, the dynamics is determined, not only by the 
orbital equations of motion for these objects, but also 
by the precessional equations for the orbital plane and 
the spin vectors themselves [6]. Note that for spinning 
compact binaries, spin effects enter the dynamics for the 
first time through spin-orbit interaction. Therefore, the 
Hamiltonian describing the leading order spin-orbit in- 
teraction can be added to the PN accurate conservative 
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Hamiltonian which gives the dynamics of nonspinning 
compact binaries. 

It is highly desirable to obtain Keplerian-type para- 
metric solution to the conservative dynamics of spin- 
ning compact binaries. This is so as 1PN accurate 
quasi-Keplerian parametrization, obtained by Damour 
and Deruelle [7], is heavily employed to analyze radio ob- 
servations of relativistic binary pulsars and to test general 
relativity in strong field regimes [8, 9]. Moreover, 2PN 
accurate Keplerian-type parametric solution for nonspin- 
ning compact binaries, obtained in Refs. [10, 11], is heav- 
ily employed to construct highly efficient "ready to use" 
search templates for nonspinning compact binaries mov- 
ing in inspiralling eccentric orbits [12]. 

Recently, 3PN accurate generalized quasi-Keplerian 
parametrization for the conservative dynamics of spin- 
ning compact binaries in eccentric orbits was obtained in 
Refs. [13, 14]. More precisely, the PN accurate paramet- 
ric solutions are available for the following two distinct 
cases: (i) the binary is composed of two compact objects 
of equal mass and two arbitrary spins, and (ii) the two 
objects have unequal masses but only one object has a 
spin. The PN accurate parametric solution may be used 
to describe the following two astrophysical situations. As 
the measured neutron star (NS) masses in NS-NS bina- 
ries usually cluster around 1.4M©, case (i) may be used 
to describe such NS-NS binaries. The case (ii) may be 
applicable in the instance of a neutron star orbiting a 
much larger, rapidly rotating black hole. In that case, 
the black-hole spin will dominate the orbital precession, 
and the neutron star spin may be neglected. 

The above parametrization, for the two cases consid- 
ered, shows that the dynamics is integrable and cannot 
be chaotic [5]. The parametric solution will be employed 
to construct "ready to use" search templates of spinning 
compact binaries in inspiralling orbits [15]. Further, it 
will also be useful to construct a fully 2PN accurate "tim- 
ing formula", required to analyze binary pulsar observa- 
tions by the Square Kilometre Array (SKA) [16]. 
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In this paper, we compute, using our parametrization, 
the observable dimensionless relativistic periastron ad- 
vance parameter k, relevant for binary pulsars. The 
expression for k that we derive will be applicable, like 
the parametrization computed in Ref. [14], only for the 
two distinct cases mentioned above. Our expression for 
k is in agreement with the one obtained by Damour 
and Schafer [10]. It should be noted that Damour and 
Schafer employed the evolution equation for the Laplace- 
Runge-Lenz vector, computed from a Hamiltonian, to 
obtain the spin-orbit contributions to the parameter k. 
In this paper, we employ our generalized quasi-Keplerian 
parametrization, derived in Ref. [14], to compute the pa- 
rameter k. 



There are two reasons that motivated us to do the 
current computation. The 2PN accurate expression for 
k, derived by Damour and Schafer, valid for arbitrary 
masses and spins, is extremely important due to the dis- 
covery of the double pulsar PSR J0737-3039 [17]. This is 
because the accurate timing of the double pulsar over a 
period of few years should lead to the measurement of the 
moment of inertia of one of the pulsars [18]. It is from the 
measurement of k to 2PN order, which includes leading 
order spin-orbit contributions, that one will be able to 
measure the moment of inertia of the fast rotating pulsar 
in this double pulsar [10]. The measurement of the mo- 
ment of inertia of a neutron star for the first time will lead 
to strong constrains on its equation of state [19, 20]. We 
compute k, as mentioned earlier, using Keplerian-type 
parametrization with the aim of presenting an alterna- 
tive way of obtaining this observable quantity. 

A second reason for computing k is the existence of 
certain intriguing features appearing in our parametriza- 
tion. We observed that the factor that modulates the true 
anomaly, due to spin-orbit interaction, contains a contri- 
bution, independent of the spins themselves, and of 1PN 
order [see Eq. (5.6k) in Ref. [14]]. We want to make sure 
that this feature will not appear in the observable param- 
eter k, computed using our generalized quasi-Keplerian 
parametrization. 

The paper is structured in the following manner. In 
Sec. II, we briefly describe the 3PN accurate Keplerian- 
type parametrization for the conservative dynamics of 
spinning compact binaries, moving in eccentric orbits, 
when spin effects are restricted to the leading order spin- 
orbit interaction, derived in Ref. [14], and point out its 
salient features. Sec. Ill and Sec. IV deal with the issues 
related to the derivation of the parameter k. We show for 
the first time how to derive k using the parametrization, 
described in Sec. II, after describing how Damour and 
Schafer obtained k. We summarize our results in Sec. V. 
Appendices A and B provide details of our computation 
and consistency checks. 



II. GENERALIZED QUASI-KEPLERIAN 
PARAMETRIZATION FOR SPINNING 
COMPACT BINARIES 

Recently, 3PN accurate generalized quasi-Keplerian 
parametrization for compact binaries in eccentric orbits 
with the first-order spin-orbit interaction was obtained in 
Refs. [13, 14] for the following two distinct cases: (i) the 
binary is composed of two compact objects of equal mass 
and two arbitrary spins, and (ii) the two objects have 
unequal masses but only one object has a spin. The un- 
derlying dynamics is fully specified by the following PN 
accurate (reduced) Hamiltonian H, written symbolically 
as 

H{r,p,S 1 ,S 2 ) = H N {r,p) + H 1PN (r,p)+H 2PN (r,p) 

+ H 3PN (r,p) + H so (r,p, S u S 2 ) , 

(2.1) 

where Hjq, -f/ipN, H 2P ^, and i?3PN are, respectively, the 
Newtonian, first, second and third PN contributions to 
the conservative dynamics of compact binaries, when the 
spin effects are neglected. The leading order spin-orbit 
coupling to the binary dynamics is given by i?so • In the 
above equation, r = 7t/(GM), r = \r\, and p = "P/^, 
where JZ and V are the relative separation vector and 
its conjugate momentum vector, respectively. The famil- 
iar symbols M and [i have the usual meaning, namely 
the total mass and the reduced mass. The explicit ex- 
pressions for PN corrections, associated with the mo- 
tion of nonspinning compact binaries, were obtained in 
Refs. [10, 21]. The spin-orbit contribution, available in 
Refs. [10, 14, 22], reads 

H S o(r,p,S 1 ,S 2 ) = ^(rxp)-S eS , (2.2) 

where the effective spin S e g is defined by 

S eff = S 1 S 1 + S 2 S 2 , (2.3) 

with 

Si = \ + \ (l - V^v) , (2-4a) 
5 2 = 1 l + \{l + ^T^i) . (2.4b) 

The finite mass ratio r\ is given by 77 = [i/M. The 
(reduced) spin vectors S\ and S 2 are related to the 
individual spins S\ and 52 by S\ = S\/{^GM) and 
S 2 = S 2 / ([iGM), respectively. We recall that 7?., T, 
Si, and S 2 are canonical variables, such that the or- 
bital variables commute with the spin variables, e.g., see 
Refs. [6, 22]. 

In Ref. [14], we obtained the Keplerian-type 
parametrization in an orbital triad (i, j, k). In this triad, 
the unit vector i is along the line of nodes, associated with 
the intersection of the orbital plane with the invariable 
plane (ex,ey), which is the plane perpendicular to the 
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total (reduced) angular momentum J, the only constant 
vector in our analysis. The unit vector fe is defined by 
k = L/L, where L = r x p is the (reduced) orbital an- 
gular momentum vector, and hence k is perpendicular to 
the orbital plane (see Fig. 1). We display below, in its 
entirety, the parametric solution to the conservative 3PN 
accurate dynamics of spinning compact binaries moving 
in eccentric orbits, in ADM-type coordinates, with spin 
effects restricted to the leading order spin-orbit interac- 
tions. The dynamical vectors of the problem are given 

by 



r(t) = r{t) cos ip(t)i(t) 

L(t) = Lk(t) , 

S(t) = Je z -Lk(t). 



r(t) am <p(t)j(t), 



(2.5) 
(2.6) 
(2.7) 



We note that in the equal-mass case (i) the (reduced) to- 
tal spin S is the sum of two arbitrary spins S i and S2 , 
and in the single-spin case (ii) S is either Si or S2. In 
our restricted cases, namely (i) and (ii), the (reduced) or- 
bital angular momentum L and the (reduced) total spin 
S precess about the conserved (reduced) total angular 
momentum vector J = L + S at the same rate [14]. This 
configuration is called simple precession in the literature, 



e.g., see Ref. [23]. 

The time dependent basic vectors (i, j, k) are explicitly 
given by 

i(t) = cos T{t)e x + sin T(t)e Y , (2.8a) 

j(t) = — cos0 sinT(t)ex + cos 6 cosT(i)ey + sin 6 ez , 

(2.8b) 

k(t) = sinG sin T(t)ex — sin cos T(t)ey + cos ez ■ 

(2.8c) 

The angle denotes the constant precession cone angle 
of L = Lk around J and is given by 



sin 



COS0 



Ssina 
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L + S cos a 
J 



(2.9a) 
(2.9b) 



where the magnitude of the total angular momentum is 
given by J = (L 2 + S 2 + 2LS cos a) 1 / 2 and a is the con- 
stant angle between L and S. 

The time evolution for r(t), (p(t), and T(t) is given in 
a parametric and PN accurate form, which reads 
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(2.10a) 


Kt . 

— — sin iv . 

c 6 


(2.10b) 


— sin 5v , 

c 6 


(2.10c) 




(2.10d) 




(2.10e) 



The PN accurate expressions for the orbital elements e t , k, and and the PN orbital functions g±t, 9m, 

fit, fat, ht, h 6 t, fitp, fatp, 94ip, g&ip, i& v , and h 6ip , in terms of E, L, S, 77, and a, are available in Ref. [14]. The mass 
dependent coupling constant Xso connects the effective spin with the total (reduced) spin via S e ff = S1S1+S2S2 = XsoS 
and is given by 

{Si = 62 = 7/8 for (i), the equal-mass case, 
<5i(r/) or 82(7]) for (ii), the single-spin case. 

We display below the 3PN accurate expression for the internal parameter k, appearing in Eq. (2.10c) [this parameter 
was denoted by k in Ref. [14], see Eqs. (5.5c) and (5.6k) there]. It consists of two parts: the first one due to the 3PN 
accurate conservative dynamics of nonspinning point masses and the second one due to the leading order spin-orbit 
interaction 



k = &3PN + k s 



(2.11) 



where 



"3PN 
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(36960 - 40000?? + 615tt 2 t] + I68O77 2 ) 



c 2 L 2 



Xso cos a 



(2.12a) 
(2.12b) 



It should be noted that k 3PN gives the observable ad- 
vance of periastron in the case of nonspinning compact 
binaries. This results from the construction of k [see 
Sec. V in Ref. [14] for details and Ref. [10]]. However, 
when spin-orbit interactions are added to the dynamics 
of compact binaries, k does not directly lead to the mea- 
surable advance of periastron, relevant for binary pul- 
sars. This parameter A: is a measure for the advance of 
periastron only inside our parametrization and has to 
be related to the observational one. Further, note the 
presence of the Xso term (which does not depend on the 
spins), appearing in fc so at the 1PN order. It is due to 
these points, namely the Xso term that appears at the 
1PN order and k does not give the observable periastron 
advance for spinning compact binaries, that we compute 
the dimensionless observable periastron advance param- 
eter k for spinning compact binaries in this paper. Note, 
that we have corrected the sign in front of Xso/3 in fc so , 
important for the proper nonspinning limit and the fol- 
lowing calculation of the observable periastron advance 
parameter k in Sec. IV. 

Before we go onto the details of computing k, let us 
list a few salient featurs of the parametrization. The pa- 
rameter e, appearing in Eq. (2.10d), is given by the New- 
tonian contribution to say e r . Note that at Newtonian 
order all three eccentricities e r , e*, and e v , are identical 
to each other and at PN orders they are related by PN 
accurate relations [see Eqs. (5.7) in Ref. [14]]. Therefore, 
it is possible to describe the binary dynamics in terms 
of one of the eccentricities. We emphasize that, while 
adapting the parametrization to describe the binary dy- 
namics, care should be taken to restrict orbital elements 
to the required PN order. We are aware that the above 
parametrization does not lead to the classic Keplerian 
parametrization simply by putting S = 0. The detailed 
explanation of the nonspinning limit and the apparent 
lack of a simple nonspinning limit is given in Sec. IV B 
in Ref. [14]. 

The above parametric solution is usually refered to as 
the genaralized quasi-Keplerian parametrization for the 
PN accurate orbital dynamics of spinning compact bina- 
ries. It extends the computations of Damour, Deruelle, 
Schafer, and Wex [7, 10, 11, 13, 24]. 

Finally, we want to emphasize that the parametric so- 
lution, given by Eqs. (2.5)— (2.12), describes the entire 
conservative PN accurate dynamics of a spinning com- 
pact binary in an eccentric orbit, when the spin effects 
are restricted to the leading order spin-orbit interaction. 



This means the parametrization consistently describes 
not only the precessional motion of the orbit inside the 
orbital plane, but also the precessional motions of the or- 
bital plane and the spins themselves. Further, the deter- 
ministic nature of the dynamics described by the Hamil- 
tonian, given by Eq. (2.1) and treated as a self-contained 
dynamical system, naturally follows from the above para- 
metric solution. 

In the next section, we develop the basics required to 
compute k from our parametrization. 



III. RELEVANT INPUTS REQUIRED TO 
COMPUTE k 

In this section, we will introduce two new triads, in 
terms of which it will be easy to obtain k. We will also 
discuss in detail how to connect these triads and hence 
the parametrization expressed in terms of these triads. 



A. The periastron-based triad for the 
parametrization 

As we are interested in computing k using our 
parametrization, it is natural that we introduce a 
periastron-based triad (a,b,k). Later, we will express 
the relative separation vector r, given by Eq. (2.5), in 
terms of this new triad. 

Let us introduce this periastron-based triad in a ped- 
agogical way. We start by writing the angle (p, given by 
Eq. (2.10c), in the following way 

<p = v + n. (3.1) 

Remember, ip is the polar angle inside the orbital plane 
measured from the line of nodes i, associated with the 
intersection of the orbital plane with the invariable plane 
(ex,ey), and v is the "true anomaly", i.e., the angle be- 
tween the (instantaneous) periastron and the position of 
the relative separation vector r. The angle IT may be 
considered as the "argument of the periastron", i.e., the 
angle between the line of nodes i and the (instantaneous) 
periastron. By comparing Eqs. (3.1) and (2.10c), we can 
easily infer the parametric evolution for II. Therefore, 
it is obvious that the angle II contains via k the secular 
contributions to the periastron advance. 

Inserting Eq. (3.1) into Eq. (2.5) and making use of 
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Figure 1: The binary geometry and definition of various angles. Our reference frame is (ex, ey, ez), where the basic vector 
ez is aligned with the fixed total angular momentum vector J = L + S , such that J = Jez- The invariable plane (ex, ey) is 
perpendicular to J. Important for the observation is the line-of-sight unit vector N, pointing from the source (compact binary) 
to the observer. We may have, by a clever choice of ex and ey, the line-of-sight unit vector N in the ey-ez plane, k = L/L is 
the unit vector in the direction of the orbital angular momentum L, which is perpendicular to the orbital plane. The constant 
inclination of the orbital plane with respect to the invariable plane is O, which is also the precession cone angle of L around 
J. The three lines of nodes are associated with the intersections of the orbital plane, the invariable plane and the plane of the 
sky. The orbital plane intersects the invariable plane at the line of nodes i, with longitude T, measured in the invariable plane 
from ex = p, which also gives the phase of the orbital plane precession. The orbital plane intersects the plane of the sky at 
the line of nodes U, with longitude fl, measured in the plane of the sky, from ex = P- The subscript i in i; results from the 
corresponding orbital inclination angle i. The dynamical longitude of periastron, measured in the orbital plane from i, is II. 
The observational longitude of periastron, measured in the orbital plane from ii, is uj. The above figure implies uo = A + II. 



familiar trigonometric relations leads to 

r = r cos(v + II) i + r sin(u + II) j 
= r cos v a + r s'mv b , (3-2) 

where the new basic vectors are given by 

a = cos II i + sin II j , (3.3a) 
b= -sinlH + cosIIj. (3.3b) 

Now, it is obvious that a is the unit vector pointing to- 
wards the (instantaneous) periastron, as r is periodic in 
the angular variable v with period 2it in these new ba- 
sic vectors a and b. Hence it is straightforward to de- 
fine the periastron-based orthonormal triad (a, b, k) = 
(a, k x a,fc), by rotating the (i,j,k) frame around the 
orbital angular momentum unit vector fc, by an angle II. 



This gives us 

a\ / cosn sinn 0\ I i\ 
b = -sinn cosn j\ , (3.4) 
k) \ 1/V fc / 

where n is the angle between i and the (instantaneous) 
periastron unit vector a, as explained above. 

The periastron-based vectors (a, b) are co-rotating 
with the precessing (instantaneous) Newtonian ellipse. 
This is similar to the behaviour of a body fixed reference 
system usually encountered in classical mechanics [25]. 
Hence, the set of angles (T,n, 0) is equivalent to the 
Eulerian angles defined there. 

In the next subsection, we will introduce an observer- 
based triad and investigate how to connect the observer 
and periastron-based triads. 



Connecting the observer and periastron-based 
triads 



The second new triad that we introduce is a fixed one 
associated with the observer on earth and denoted by 
(p, q, N). In this fixed triad p and q are two orthogonal 
unit vectors in the plane of the sky, i.e., the plane trans- 
verse to the radial direction TV, pointing from the source 
(compact binary) to the observer (see Fig. 1). 

In order to compute k using our parametrization, it is 
natural to expect that we need to connect the periastron- 
based triad (a, b, k) to the observer triad (p, q, N). It is 
possible to connect these two frames with the help of ro- 
tational matrices quite similar to the way in which the 
(i,j,k) frame was connected to the (p,q,N) frame in 
Ref. [14] [see Sec. VI in Ref. [14] for details]. In this 
paper, we come from the observer triad (p, q, TV) to the 
periastron-based triad (a, b, k) with the help of four suc- 
cessive rotations: first, a rotation around p via the con- 
stant angle iq to come to the reference triad (ex,ey, ez), 
second, a rotation around ez by an angle T, third, a rota- 
tion around the direction of the line of nodes i [associated 
with the intersection of the orbital plane with the invari- 
able plane (ex,ey)] by an angle 6, and, finally fourth, 
a rotation around the orbital angular momentum direc- 
tion k by an angle II. In terms of rotational matrices — 
denoted by 7?.[rotation angle, rotation axis] — we have 

'o(f)\ 

b(t) ^n[u(t),k(t)}n[e,i{t)}TZ[T{t),ez} 
K k(t)J 



p 

x K[io,p] | Q 



(3.5) 



where the first rotation appears at the right-hand side, 
the second one appears left from the one before, and so 
on. 

We note, that the approach to generate the (a, 6, k) 
frame via these rotations may be also adapted for scenar- 
ios, where the angle depends on time. The constant 
behaviour of is a direct consequence of choosing two 
distinct cases (i) and (ii), where our parametrization is 
valid and hence we will use the fact that is constant 
only at the end of our analysis. 

It should be noted that there are other possible ways to 
connect the orthonormal triads (p, q, TV) and (a, 6, k). A 
second way to generate the (a, b, k) frame is described in 
great detail by Damour and Schafer in Sec. 5 of Ref. [10] 
and involves three rotational matrices. This way is closely 
linked to the observation as it involves rotational matri- 
ces defined with the help of the observational longitude of 
periastron uj and the orbital inclination angle i. Damour 
and Schafer connected the observer triad (p, q, TV) to 
the periastron-based triad (a, 6, k) via the following three 
successive rotations: first, a rotation around TV by an an- 
gle f2, second, a rotation around the direction of the line 
of nodes ii [associated with the intersection of the orbital 




plane with the plane of the sky (p, q)] by an angle i, and, 
finally third, a rotation around the orbital angular mo- 
mentum direction k by an angle u>. The corresponding 
rotational matrices read 

= U[uj(t),k(t)} K[i(t),ii(t)] K[fl(t),N] I q 

(3.6) 

We note, that for the unit vector ii along the line 
of nodes, associated with the intersection of the orbital 
plane with the plane of the sky, the subscript i is em- 
ployed as ii is related to the orbital inclination angle i 
and to avoid any confusion with the unit vector i. 

In the next section, we obtain k with the help of the 
parametrization and the ideas discussed in this section. 



IV. THE ADVANCE OF PERIASTRON 
RELEVANT FOR BINARY PULSARS 

In this section, we will compute fc, the observable di- 
mensionless periastron advance parameter due to k. Fol- 
lowing Damour and Schafer, the PN accurate total ob- 
servable dimensionless periastron advance parameter k 
can be written as 



^3PN 



(4.1) 



where fc 3 pN is the usual PN accurate periastron advance 
parameter up to 3PN order, corresponding to the non- 
spinning case [see Refs. [10, 13] for the 2PN and 3PN 
contributions, respectively], and fc so denotes the contri- 
bution due to spin-orbit interaction. It should be kept in 
mind that the 3PN contribution to /«3pn will not be re- 
quired while probing the dynamics of relativistic binary 
pulsars even with SKA. 

Let us first briefly summarize how Damour and Schafer 
obtained k, especially k so , in the next subsection. 



A. The derivation of fc so — The approach of 
Damour and Schafer 

Timing of binary pulsars allows via the observed "Roe- 
mer time delay" the direct measurement of the argument 
of periastron uj [8]. Therefore, Damour and Schafer de- 
fined the relevant secular periastron advance as (d,Lu/dt) t , 
or equivalently the timing parameter (not to be confused 
with the orbital angular-momentum unit vector fe) 



1 /duj\ 
~dt/ t ' 



(4.2) 



where n is the mean orbital motion that appears in the 
PN accurate Kepler equation, given by Eq. (2.10b), and 
(. . .) t denotes the orbital average. 

Damour and Schafer computed the contribution to 
diu/dt due to spin-orbit interaction in the following way. 
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First, they computed the precessional equations for the 
orbital angular momentum vector C and the Laplace- 
Runge-Lenz vector A., which tracks the periastron ad- 
vance, with the help of the 2PN accurate version of 
Eq. (2.1). These precessional equations where orbital av- 
eraged to obtain the rotational velocity vector, denoted 
by ft in Ref. [10]. The expression for doj/dt was extracted 
from ft by decomposing ft along three non- orthogonal 
directions [see Eq. (5.11) in Ref. [10]]. In this way a 
parametric solution to the dynamics of spinning compact 
binaries was not required to compute the 2PN accurate 
measurable dimensionless periastron advance parameter 
k in Ref. [10], and hence the formula of Damour and 
Schafer for k is more general. 

In the next subsection, we will compute dto/dt us- 
ing our Keplerian-type parametric solution and the two 
different ways of connecting the (p, q, N) frame to the 
(a, 6, k) frame. 



B. The derivation of fc so — 
approach 



The parametric 



This subsection deals with the computation of doj/dt, 
arising due to spin-orbit interaction, using our Keplerian- 
type parametrization. 

Following Damour and Schafer, we introduce the ro- 
tation velocity vector ft, containing all involved angular 
velocities, and extract dto/dt from it in the following way. 
We have already noted that there are two different ways 
to connect the observer-based triad to the pariastron- 
based one [see Eqs. (3.5) and (3.6) and related discussions 
there]. The first way, given by Eq. (3.5), involves the an- 
gular variables 11(f), 6, T(f), and io, appearing in our 
parametrization, and the second way, given by Eq. (3.6), 



involves the angles u)(t), i(t), and ft(t), where two of them 
are relevant for the observations, namely u(t) and i(t). 
Therefore, we can decompose the rotation velocity vector 
n in the following two ways: 



dll, 
——k - 

dt 


de . 


dT 




' ~dt 


duj , 


di . 


dfi 




dt 



dip 
dt 1 



(4.3a) 
(4.3b) 



Note, that the above two different decompositions 
are along four/three non- orthogonal directions, respec- 
tively. The elementary rotation velocities, appearing 
in Eq. (4.3a), can be defined with the help of our 
parametrization, whereas in Eq. (4.3b) we have doj/dt: 
the rate of change of the argument of the periastron. 
Further note, that Eq. (4.3b) is identical to Eq. (5.11) of 
Ref. [10]. 

In order to extract duj / dt and to relate it to the angular 
variables associated with our parametrization we proceed 
as follows. Since all vectors, appearing in Eqs. (4.3), can 
be expressed in terms of (a,b,k), it is clear that the 
rotation velocity vector ft may also be written as the 
sum of three elementary rotation velocities, along three 
orthogonal directions a, b, and fe, namely 



ft 



«(o)« 



(4.4) 



It is then straightforward to compute the components 
^(o)j ^(6) j an d fi(fc) m terms of the elementary rotation 
velocities and angles appearing in Eqs. (4.3). The com- 
ponents of fl in the co-rotating periastron-based (a, b, k) 
frame — using the parametrization-based variables and 
the observation-based variables, respectively — are given 
by 



de dT 
"fal = ~ t cos II H — — sin 6 smll 

w dt dt 



—r- (cos II cos T — cos 9 sin II sin T) = — coscj 
dt y ' dt 



dfl 
~dl 



sm i sin oj . 



de dr 

SZfhi = — sin II H — — sin cos II 

w dt dt 



—r-(— sin II cos T — cos 9 cos II sin T) = — - sinw 
dt y 1 dt 



dfl 
~dt 



sin i cos w . 



fi 



(fe) 



dH dT dig ~ „ du> dfl 

— — I — — cos 9 H — — sm 9 sm T = — — | — — cos i . 
dt dt dt dt dt 



r 



(4.5a) 
(4.5b) 
(4.5c) 



These general set of equations provides the transforma- 
tion from the parametrized angles to the angles of interest 
for the observation. We note, that the last component, 
Cl(k) , gives the so-called intrinsic contribution to the pe- 
riastron advance t ^ mtrmsic ) defined by w lntrinsic = ft k. 
This is the component of ft acting within the orbital 
plane. However, in the timing of binary pulsars the pa- 
rameter k enters in the (relativistic) argument of the pe- 
riastron u>, and so we have to extract dwjdt from the 
above set of equations. 



A close inspection of Fig. 1 indicates that io is constant, 
since it connects two invariant reference frames. For the 
two distinct cases (i) and (ii) , where our parametrization 
is valid, we also have 9 as a constant angle. This allows 
us to come from Eqs. (4.5) to the reduced set of equations 
for the components £l( a ), ft(u\, and fl(k), which reads 



dT 

~dt 
dT 

~dt 



sin 9 sin II 
sin 9 cos II 



di dfl . . . 

— cos u) -\ — — sm i sm u . 
dt dt 



di 
"dt 



dfl 
~dt 



sm i cos oj . 



(4.6a) 
(4.6b) 
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dU dT duj dn 

— — I — — cos = — — I — — cos i . 
dt dt dt dt 



(4.6c) 



Solving the last equation, namely Eq. (4.6c), for dui/dt 
leads to 

dw dU dT „ dtt , d N 

— = — + — cos 9 — cosi. (4.7) 

dt dt dt dt y J 

The first and the second contribution to dw/dt are di- 
rectly provided by our parametrization. We also figured 
out that they can be obtained via the time evolution of 
the Laplace- Runge-Lenz vector without any parametriza- 
tion (see Appendix B for details). The third contribution 
to dbj/dt in Eq. (4.7), namely (dCt/dt) cost, needs a closer 
inspection. To obtain d&l/dt, we multiply Eq. (4.6a) with 
s'mu! and add it to the product of Eq. (4.6b) times cos u>, 
which results in 



dfl 



1 dT 



dt sini dt 



sinG (sin II sin w + cos II cos a; ) . (4i 



Using now the explicit expressions for sinw and cosu, 
given by Eqs. (A9) in Appendix A, in the above equation, 
we get 



dn 



1 dT 



dt sin 2 i dt 



Se(C io S@ + S io C @ cosT), (4.9) 



where Sq , Cq , £j , and Cj are shorthand notations for 
the constant trigonometric functions sin0, cosG, sinzo, 
and cosio, respectively. The time evolution for the or- 
bital inclination i, appearing in Eqs. (4.7)-(4.9), is simply 
given by Eqs. (A4) in Appendix A. 

Let us now explicitly compute dH/dt, dT /dt and hence 
dCl/dt with the help of our parametric solution. Using 
Eqs. (2.10c), (2.11), and (3.1) for II, and Eq. (2.10d) for 
T, we display below the relevant parametric expressions 
symbolically 



II - n = f&3PN + fcso) v + (■•■) sin 2v + 
+ (...) sin4v +(...) sin5v , 



.) sin 3v 
(4.10a) 
(4.10b) 



The structure of Eqs. (4.10) indicates that we can handle 
the PN and the spin-orbit induced contributions to dw/dt 
separately. 

To obtain the spin-orbit contributions to {dui/dt} t , we 
compute with the help of Eqs. (4.10) and (2.12b) the or- 
bital averaged spin-orbit (so) contributions to the above 
required time derivatives 



dn\ sc 

dt/t 



i / 1 



dT 
~dt 



l/i 



X BO (-i - ZScosu) , (4.11a) 



X so (i + 5cosa), (4.11b) 



'dfi \ 1 / 1 \ r, . cosi 



x (C l0 S e + 5 l0 C e cosT) 



(4.11c) 



The quantity (l/r 3 ) t gives the orbital average of l/r 3 , 
r being the scaled radial separation, to the Newtonian 
order [see, e.g., Ref. [10]] 



1 



1 



3 / t a 3 ( 1 _ e 2)3/2 (l_ e 2)3/2 



(4.12) 



With the help of Eqs. (4.7), (4.11), and (4.12), it is 
straightforward to obtain the spin-orbit induced contri- 
bution to k, defined by Eqs. (4.1) and (4.2), and it reads 



D 



-I cos a — sin a : „ (Gj oq 
sin i 



S ia Ce cos T) 



(4.13) 



where B is given in terms of our reduced variables by 

nXsoS 



B = 



c 2 (l - e 2 ) 3 / 2 



(4.14) 



The proper nonspinning limit of Eqs. (4.13) and (4.14), 
can be easily deduced as when S — > 0, J is identical 
to L, which implies 6-^0 and io — ► i, i being the 
angle between N and L, and therefore k so vanishes. The 
disappearance of fc so as S — * indicates that the presence 
of the 1PN accurate x so term in fc so , given by Eq. (2.12b), 
does not have any observational consequence. 

Our result for fc so is in agreement with Eqs. (5.17) and 
(5.17a) of Ref. [10], when we specialize these equations 
to the simple precession (for the equal-mass case and for 
the single-spin case) 



B 



sin 2 i 



[(1 -3sin 2 i)(fc-s) -cosi(N-s)] , (4.15) 



where s = S/S, k ■ s = cos a and the coefficient B is 
given here in terms of the non-reduced variables, used by 
Damour and Schafer in Ref. [10] and in Appendix B, by 



B = 



c 2 M(l - e 2 ) 3 / 2 



(4.16) 



The dimensionless mass dependent coupling constant A so 
is related to x so by A so = Xso/v (see Appendix B for 
details). To check the above agreement between our re- 
sult and the one of Damour and Schafer, we have used 
Eqs. (2.9) and Eq. (A4a). The two expressions for the 
coefficient B, namely Eqs. (4.14) and (4.16), are identi- 
cal. 

Let us now briefly consider the nonspinning 3PN ac- 
curate contributions to k. These contributions, enter 
only via II, orbital averaged analogue to Eq. (4.11a), our 
computation, and result in &3pn = However, it 

should be noted that the expression for fc.3PN, as given 
by Eq. (2.12a), is not yet ready for application. This is 
so as this expression for fc 3PN depends on the conserved 
theoretical quantities E and L rather than on observable 
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quantities like mean motion n and proper-time eccentric- 
ity &r- However, we have checked that up to 2PN order 
our expression for /c 3 pn is in agreement with the one of 
Damour and Schafer, obtained in Ref. [10]. 

Therefore, the 3PN accurate observable periastron ad- 
vance parameter k becomes 



^3PN 



(4.17) 



where /c 3 pn is the orbital 3PN contribution, correspond- 
ing to the nonspinning case, given by the 3PN accurate 
contribution to Eq. (2.11), namely Eq. (2.12a). The spin- 
orbit contribution k so is given by Eqs. (4.13) and (4.14) 
or equivalently by Eqs. (4.15) and (4.16). We recapitu- 
late that fc 3 pN = fc 3 pN and fc so ^ k so , where fc 3 pN and 
fc so are the corresponding contributions to k, as given 
by Eqs. (2.11) and (2.12). This results from the con- 
struction of k in Ref. [14] by combining the parameter 
^3PN = ^3PN of the nonspinning case with the internal 
parameter fc so , given by Eq. (4.45e) in Ref. [14], of the 
spin-orbit parametrization [see Sec. V in Ref. [14] for de- 
tails] . 

Finally, let us present the following geometric explana- 
tion why k k. The difference between k and k results 
from the fact, that the angle II — associated with the pa- 
rameter k — is measured from i and in contrast to that, 
the angle u> — associated with the observable parameter 
k — is measured with respect to ii. The difference be- 
tween u> and II is denoted by A in Fig. 1 (see Appendix A 
for details). In addition, we want to point out, that the 
last rotation TZ[u>(t), k(t)] in Eq. (3.6) can be decomposed 
into two parts 

Tl[u{t), k(t)} = K[IL(t), k(t)] K[A(t), k(t)} , (4.18) 

and so it is obvious that w(t) = A(t) + H(t) and corre- 
spondingly we have duj/dt = dA/dt + dll/dt. The con- 
tribution of dA/dt causes the difference between k so and 



V. CONCLUSIONS 

We have presented a new method to compute the PN 
accurate observable dimensionless relativistic periastron 
advance parameter A: of a spinning compact binary us- 
ing our Keplerian-type parametric solution to the PN 
accurate conservative dynamics of spinning compact bi- 
naries moving in eccentric orbits [14]. The expression for 
k is applicable, like our parametric solution, only for the 
following two distinct cases, namely, case (i), the equal- 
mass case and case (ii), the single-spin case. Our result 
is consistent with a more general formula for k, derived 
by Damour and Schafer, using a different procedure [10]. 
The computation also helped us to clarify certain intrigu- 
ing features that we observed in our parametrization. 
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Appendix A: SOME ASPECTS OF THE BINARY 
GEOMETRY AND RELATED DYNAMICS 



In this section, we make contact with the variables and 
angles, commonly used in the literature [10, 26]. It is ob- 
vious that we have three planes in the binary geometry, 
namely the orbital plane, the invariable plane and the 
plane of the sky. The associated orthonormal triads are 
(i,j,k) } (ex,ey,e^), and (p, q, N). Therefore, we may 
define three lines of nodes and associated inclination an- 
gles in the following way. We infer from Fig. 1 that the 
plane of the sky meets the invariable plane along a line 
of nodes, characterized by p = e x , as we choose the line- 
of-sight unit vector TV to lie in the ey-ez plane. The 
associated constant angle io between N and J = Jez 
may be treated as the inclination angle, related to the 
above line of nodes. In the second case, we consider the 
line of nodes i, created by the intersection of the or- 
bital plane with the invariable plane. The corresponding 
inclination angle, namely the angle between the orbital 
angular momentum unit vector k = L/L and J = Jez, 
is denoted by 6. Finally, we consider the intersection of 
the orbital plane with the plane of the sky. We associate 
ii with this line of nodes and the related inclination angle 
i, the angle between k and the line of sight unit vector 
N, is the usually defined orbital inclination i and hence 
the use of the subscript i in ii. 

It is easy to figure out that the unit vectors, associated 
with the last two lines of nodes, i and ii are defined by 



ii = 



ez x k 
\ez x k\ 
N xk 
IN x fcj 



(Ala) 
(Alb) 



Using fe, expressed in terms of {ex,ey ,ez) and 
(P,Q,N), given by 

k = Sq sin Tex — Se cos Tey + C^ez 
= SesinTp-|- (-S io C @ - C io S@ cosT)q 
+ (C io C e -S io S e cosT)JV, (A2) 



we compute the following expressions for i and i, as 

(A3a) 



i = cos Tex + sin Tey 
= cos Tp + d sin Tq + Si sin TN . 
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i% = -. — r [(Si Ce + d Se cos T)p + Se sin Tq] , 

(A3b) 

where sini = |7V x k\. The time evolution for the (ob- 
servational) orbital inclination i is simply given by 



cos i = N ■ k = Ci Cs — Si Se cos T , 
sinz = \N x k\ 

= [(S io C @ + C l0 S e cos T) 2 + sin 2 T] 



[I -(C io Ce-S io S e cos?) 



21 1/2 



(A4a) 



1/2 



(A4b) 



The equation for cosi is in agreement with Eq. (6) in 
Ref. [26]. We may also obtain a differential equation for 
the time evolution of i, as given in Ref. [10]. To get 
that, we differentiate cosi = N ■ k with respect to time 
and make use of k = L/L with L = S e ff x L/(c 2 r 3 ), 
sini = |7V x fe| and Eq. (Alb). In this way we deduce the 
following differential equation for the orbital inclination 



di 



1 N x k 

rS e ff ' ; — — 



(A5) 



which agrees with Eq. (5.15) of Ref. [10], whereas we 
note that the quantity characterizing the ascending node 
in Ref. [10], namely i, is identical to our ii. 

Following Ref. [10], let fl be the longitude of the line 
of nodes z, , as measured from p in the p-q plane (plane 
of the sky), which implies 



ii = cos flp + sin fl q . 



(A6) 



The time evolution for fl is defined by either of the fol- 
lowing two equations 

cos Q = p ■ ii = -r—(Si a C® + Ci a Se cos T) , (A7a) 
smi 



sin fi = \p x ii | = — — Se sin T . 

sin i 



(A7b) 



Finally, we note that it is also possible to obtain the 
time evolution for the longitude of periastron lu, mea- 
sured from ii in the orbital plane, as defined in Ref. [10]. 
To do that, it is advantageous to define a periastron- 
based orthonormal triad (a, b, k) = (a, k X a, k). This is 
easily obtained by a rotation of the (i,j, k) frame around 
the orbital angular momentum unit vector fc, by an angle 
II. This gives us 



-(C io S e + S^Ce cos T) cos II] 
1 



(A9a) 



sinw = \ii x a\ = — — : [(d Se + Si Ce cos T) sin IT 



-Si sin T cos II] 



(A9b) 



The parametric solution to II, which also gives its time 
evolution, is given by our parametrization via Eq. (3.1). 
In addition, the parametric solution to II can be also 
obtained by noting that a = A./\A\, where A- is the 
Laplace-Runge-Lenz vector given by 



A = V x (TL x V) - ^ 2 GM^ 



= f i 2 GM 



p x (r x p) 

r 



(A10) 



Further details are described in Appendix B. 

Eq. (A9b) for sincj is in agreement with Eq. (6) in 
Ref. [26]. The additional check sin 2 lu + cos 2 lu = 1 leads 
to 



sini = [(C io S e + S io C e cosT) 2 + S 2 sin 2 T] 



1/2 



(All) 



which is equivalent to Eq. (A4b) (using trigonometric 
relations) . 

Since ii, i, and a lie in the orbital plane, we may define 
A such that it is the angle between i t and i. This implies 
that lu = II + A (see Fig. 1) and the time evolution for 
A is given by 



cos A 
sin A = \ii x i 



ii-i = (C io S e + S io C e cos T) 

SIM 



-S in sinT . 



(Al2a) 
(Al2b) 



Finally, let tp-os be the polar angle of r, measured from 
ii, as defined by Damour and Schafer in Ref. [10]. As we 
define (p from our i, this leads to <^ds = </> + A. Using 
these inputs, it is possible to write the relative separation 
vector r in (a, b) in terms of <£ds and lu as 

r = r cos((fBS — w)a + r sin((^DS — w )b , 
= rcosva + rs'mvb, (A13) 

This is consistent with Eqs. (5.5b), (5.6), and (5.10) in 
Ref. [10]. In our angular variables, Eq. (A13) reads 

r = r cos(</j — II)a + r sin((^ — II)b , 
= rcosva + r sinvb. (A14) 



a\ / cosll sinll 
b = -sinn cosIT 
k \ 1 




(A8) 



where II is the angle between i and the (instantaneous) 
periastron unit vector a. The time evolution of lu is given 
by either of the following relations 



cos lu 



[—Si sin T sinll 



Appendix B: THE COMPUTATION OF T AND n 
VIA THE LAPLACE-RUNGE-LENZ VECTOR 



In this section, we present an elegant way, starting 
from the famous Laplace-Runge-Lenz vector, to obtain 
the time evolution of the apsidal motion without any 
parametrization. The quantities in this section are non- 
reduced for a simple comparism with the work of Damour 
and Schafer [10], and the work of Wex and Kopeikin [27]. 
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We define the following functions of the canonically 
conjugate phase space variables TZ and V 



€(71, V) = KxT. 



(Bla) 



TZ 



A(K, T) = V x (H x V) - fSGMjj- . (Bib) 

The first equation gives the orbital angular momentum 
vector C and the second one the apsidal vector A, which 
is the famous supplementary first integral of the Kepler 
problem. The vector A is usually named after Runge 
and Lenz or Laplace, though he was first discovered by 
Lagrange. 

These vectors evolve according to the fundamental 
equations of Hamiltonian dynamics 



dC 
~dT 
dA 
~dT 



{C,H}, 
{A, H} , 



(B2a) 
(B2b) 



where {...,...} denotes the Poisson brackets and H is 
the full Hamiltonian for the relative motion of our bi- 
nary system. As C and A are first integrals of Hn, only 
TtiPN + ^2PN + ^3PN + Ti-so contribute to the right-hand 
sides of Eqs. (B2). The effects of 7Yipn + W 2 pn + H3PN 
are already studied in the literature, and so we just want 
to consider the contribution of 



n so (n,v,s u s 2 ) 



G 

^3 



c s 



eff • 



where the effective spin <S e ff is given by 



<5 e ff = I 2 



3m 2 \ 



3mi 
2m 2 



(B3) 



(B4) 



We recall that Tt, "P, Si, and <S 2 are canonical variables, 
such that the orbital variables commute with the spin 
variables, e.g., see Refs. [6, 22]. Using the basic Poisson 
brackets {K 4 ,Vj} 
s ijk S$, jyc.vc.'! = 

{R\Si} = {Vi,S{} = {Vi,Si) = 0, we find 



8), {Si,S(} = e^S*, {S l 2 ,Si} 



{£,Hso} = 
{AH so } = 



G 

c 2 n 3 

G 
c?-R? 



S eS x C , 



(B5a) 



Tt 

S ee x A-3(C ■ S eS )C x — 



(B5b) 



After orbital averaging these expressions, we can write 
down the equations for the precessional motion, induced 
by the spin-orbit interaction, in the vectorial form with 
the help of the rotational velocity vector fl so 



dT I T 



xC. 



(B6a) 



dAV° = 
dT/, : 



— l~2 so x A , 



where 



— /— \ 



»efF 



(B6b) 



(B7) 



and {. . .)t denotes the orbital average. 

In case of the so-called simple precession, where the 
orbital angular momentum C and the total spin S pre- 
cess about the conserved total angular momentum vector 
J = C + S at the same rate [14, 23, 27], the precession of 
the orbital plane and the longitude of periastron are best 
described by the angles T and II, defined with respect 
to the invariable plane, i.e., the plane perpendicular to 
J = Jez (see Fig. 1 for details): 



dC dT 

dT = dT ezXC ' 



dA 

dT 



dT 
dT 



ez 



dn 
dr' 



k x A, 



(B8a) 
(B8b) 



where k = C/C denotes the unit vector in the direction of 
the orbital angular momentum vector C The conditions 
for the simple precession are fulfilled in our restricted 
cases, namely (i), the equal-mass case and (ii), the single- 
spin case. In these cases the effective spin «S e ff and the 
total spin S are related by <S e ff = A S0 <S, where the mass 
dependent coupling constant is given by 



2 + 3/2 = 7/2 
2 + 3m 2 /(2mi) 
2 + 3mi/(2m 2 ) 



for (i), the equal-mass case, 
for (ii),Si^0,5 2 = 0, 
for (ii),Si = 0,5 2 ^ 0. 



Note, that A so = Xso/V- Using S e ff = A S0 5 in Eq. (B7) 
and comparing Eqs. (B6) with an averaged version of 
Eqs. (B8), we obtain the precession in terms of the angles 
T and II 



dT\ s ° _ G 
dT/, ~ ~- 



rp O ^ V\L J rj-t 

mv° _ g 



dT / , 



TV I, 



A s „£ -1-3 



CS 



(B9a) 
(B9b) 



We note, that in case of dT /dT it is even possible to 
obtain a non-averaged equation, if one uses Eq. (B5a) 
instead of Eq. (B6a) for the comparism with Eq. (B8a). 

Our results, given by Eqs. (B9), are in agreement with 
Eqs. (4.11) and with Eqs. (38) and (42) in Ref. [27]. 
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